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The gravitational collapse of a triplet scalar field is examined assuming a hedgehog ansatz for the 
scalar field. Whereas the seminal work by Choptuik with a single, strictly spherically symmetric 
scalar field found a discretely self-similar (DSS) solution at criticality with echoing period A = 3.44, 
here a new DSS solution is found with period A = 0.46. This new critical solution is also observed 
in the presence of a symmetry breaking potential as well as within a global monopole. The triplet 
scalar field model contains Choptuik's original model in a certain region of parameter space, and 
hence his original DSS solution is also a solution. However, the choice of a hedgehog ansatz appears 
to exclude the original DSS. 
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I. INTRODUCTION 

Quite an industry has developed around the study of 
nonlinear phenomena occurring at the threshold of black 
hole formation first discovered by Choptuik Q| . A num- 
ber of reviews detail the substantial progress in under- 
standing these phenomena [2|,|j. By specifying a form 
of initial data parameterized by some parameter p, nu- 
merical evolutions can determine whether the value of 
p is sufficient to form a black hole. Such a solution is 
considered super-critical, while for small p, no black hole 
forms and the solution is sub-critical. Separating super- 
critical solutions from sub-critical solutions is the critical 
solution parameterized by the critical value p* . In anal- 
ogy with statistical mechanics where the black hole mass 
represents the order parameter, cases where infinitesimal 
mass black holes are formed are called Type II, and cases 
in which only finite mass black holes are formed are called 
Type I. 

Here, I present a new Type II critical solution occur- 
ring within a model with a triplet scalar field. By as- 
suming a hedgehog ansatz for the triplet, I find a crit- 
ical solution which is, in a strict sense, non-spherically 
symmetric. Adding a potential allows the study of crit- 
ical phenomena within a global monopole. No Type I 
critical behavior around the static monopole solutions is 
observed. 

As a generalization of Choptuik's original model, the 
triplet scalar field contains the single scalar field model 
within a certain region of parameter space. However be- 
cause of the different assumed symmetries of the fields, 
the hedgehog ansatz excludes this region. Hence, it is not 
surprising that I find a new solution, one which is dis- 
cretely self-similar (DSS) with echoing period A = 0.46, 
much smaller than that for either the original DSS (A — 
3.44) or that of the Einstein Yang-Mills (EYM) model 
(A = 0.74) |. 

An ansatz similar to this one (an SU{2) valued 
hedgehog ansatz) was used to study the collapse of 
Skyrmions ^-Q]. Both these models allow the con- 



sideration of non-spherically symmetric initial data, al- 
beit quite special data, within the comparatively simple 
framework of one spatial dimension, namely a radius r. 

In Section ||, I present the most general form of the 
model under consideration here. First, I choose a triplet 
scalar field accompanied by the usual symmetry- 
breaking potential with vacuum value |$°| = 77 and cou- 
pling A. For certain evolutions discussed later, I also in- 
troduce a free, massless scalar field ip{r,t) which couples 
to only through gravity. 

With the most general form of the equations of motion. 
Section HI then discusses various aspects of reg ions of the 
parameter space. In particular. Section III A| presents a 
new critical solution obtained with the hedgehog ansatz, 
and Section IIIB discusses the stability of this solution 
with respect to the addition of ip, a free scalar field. Sec- 
tion III C| addresses the results obtained in the interior of 
a monopole. I then conclude in Section IV. 



II. THE MODEL 

Letting represent an S'0(3)-valued triplet scalar 
field, and ip a single-valued, free, massless scalar field, 
the Lagrangian for the model is 



m-ri^ , (1) 



where A is the coupling to the symmetry-breaking po- 
tential, and Tj is the scale of symmetry-breaking. The 
stress-energy tensor is then 



9^ U--,p^-.^ + lx 



and the equations of motion for the matter fields are 



(2) 



Dijj = 

= A*" ($°) 



(3) 
(4) 
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The hedgehog ansatz is assumed for the triplet scalar 
field 



sin 6 cos ip 

= /(r, t) f = /(r, t) ( sin^sini^ 
cos 9 



(5) 



such that the triplet field is defined to be a vector in in- 
ternal space of magnitude /(r, t) pointing everywhere in 
the radial direction, f. None of the fields <I>° is a func- 
tion only of r and t, but the magnitude / is. Thus, while 
highly symmetric, the model is not spherically symmet- 
ric in the strictest sense. With this hedgehog ansatz, 
a monopolc of unit charge is obtained by enforcing the 
boundary condition /(cxD,t) — rj. When f{oo,t) = 0, no 
monopole charge exists. 

I choose the spherically symmetric metric 



ds^ = -a^{r, t) dt^ + a2(r, t) dr^ + 
and introduce auxiliary variables 



n.0 = — 0, 
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where an overdot denotes d/dt and a prime denotes d/dr. 
The full system of equations then becomes 
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In addition to these equations of motion, appropriate 
boundary conditions are needed. At the origin of co- 
ordinates, regularity is enforced with the conditions 



$^(0,t) = 
/(0,t) = 
a(0,t) = 1 
a'(0,t) = 0. 



1>/(0,0 = 
a'(0,t) = 



(16) 



Further, a mass aspect ratio m{r, t) is defined by associ- 
ating the grr metric components of (0) and Schwarzschild 



= 1 - 



2m 



(17) 



The quantity m(r, t) then describes the amount of mass 
contained within the radius r at time t. 

At large radius, an outgoing radiation condition is ap- 
plied to the scalar field tp. For non- vanishing potential, 
the field f{r,t) has an effective mass which makes find- 
ing a good outer boundary condition difficult. While the 
use of an outgoing radiation boundary condition worked 
well in some situations, in general the value of / at the 
outer boundary was simply held constant which sufficed 
whether or not a monopole was present. Refiection off 
this boundary was often significant so that the maxi- 
mum radius was kept relatively large. This prevented 
reflections from affecting the dynamics near the origin. 



III. RESULTS 

The field equations (|9|-|l5|) are then integrated with a 
second-order accurate, Crank-Nicholson finite difference 
code. The generated solutions are found to be second- 
order convergent, stable, and mass conserving. Because 
either A or can be rescaled away, I choose either A = 0.1 
or A = for all that follows. As almost all other studies 
of Type II critical behavior have, I have made use of 
an adaptive mesh. In particular, I incorporated finite 
difference approximations to the equations (^-15) into 
the adaptive infrastructure developed by Choptuik 



A. A New DSS Critical Solution 

As the Lagrangian presented in (^ is quite general, I 
first discuss results obtained with vanishing scalar field 
ip. Specifically, if ip{r,t) is initialized initially to zero 
with vanishing time derivative it will remain vanishing. 
This simplification allows the study of the collapse of the 
triplet field alone. 

I also postpone discussion of collapse within a 
monopole until Section [II C. Hence, in this section no 
monopole is present and /(oo, t) = 0. 

I find as the threshold of black hole formation is ap- 
proached, a discretely self-similar solution is found with 
period A — 0.46. A solution which approaches this pre- 
cisely critical solution is shown in Fig. |l|. In the limit of 
precise criticality, the solution is expected to exhibit an 
infinite number of echoes and to be singular at the origin 
at the time of collapse, occurring at a finite proper time. 

The discrete nature of the self-similarity is demon- 
strated in Fig. H. This solution also exhibits the usual 
mass scaling relationship in the super-critical regime, 
demonstrated in Fig. pi 



2 







1 1 1 1 1 1 1 


1 1 1 1 1 L 


0.1 








0.05 




f \ 















4- 1 


1 1 1 V_--r 1 1 1 











1 1 1 1 1 














ll 





0.06 
0.04 
0.02 


0.08 
0.06 

0.04 
0.02 




ln(r) 

FIG. 1. Marginally sub-critical solution at late time for 
= 0. The top frame displays the hedgehog profile /, the 
middle frame shows the quantity dm/dr explained in the text, 
and the last frame shows the value of 2m/r. 
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FIG. 2. Demonstration of scale-periodicity of the new 
critical solution with A = 0.46. Here, the field / is shown at 
four different times for a marginally sub-critical solution. The 
times and the value of A were chosen to minimize the differ- 
ence /(ln(r), To) — /(ln(r) -I- A, Ti). Independently, the value 
of A was also computed by A = ln((T* - r„)/(T* - r„+i)) 
where T represents the central proper time and T* the esti- 
mated critical time of collapse. The other two profiles are then 
shown shifted in ln(r) by nA. The agreement demonstrates 
the scale periodicity. 
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ln((p-p*)/p") 
FIG. 3. Demonstration of black hole mass scaling for 
ry = 0. The circles indicate the mass of black holes formed ver- 
sus a normalized distance from criticality. The line indicates 
a least-squares fit for slope 7 = 0.119 where A^bh oc \p — p*\''. 

Numerical results indicate that variations of ry in the 
interval [0, 0.5] do not alter the nature of the critical so- 
lution. In other words, the discrete nature remains with 
A = 0.46. The fact that the triplet field's value for large 
radius vanishes while the true vacuum lies at 77 implies 
that the mass of the spacetime is divergent with increas- 
ing r. Hence, the profiles of m' and 2m jr as shown in 
Fig. |l| are not the same as when 77 = 0, differing for large 
r. However, the phenomena associated with approaching 
criticality occur for ever decreasing r, and in this regime, 
the solution remains discretely self-similar with period 
A — 0.46. This result is consistent with the potential 
not playing a role in the kinetic-energy dominated critical 
behavior. In the language of Q, the potential becomes 
"asymptotically irrelevant." 

This divergence of the mass is not as problematic as 
it might initially appear because the numerical grid is 
artificially limited to some finite value of r. The mass is 
therefore limited, but this cut-off should not affect results 
for large r because of the distance needed for any effect to 
travel from the outer boundary. By comparing otherwise 
identical evolutions with different maximum radius, no 
such effect was seen. 

It is interesting to note however that when 7^ and 
/(oo, t) = 0, the value of 2m/ r was seen to be divergent. 
This situation represents an infinite mass spacetime (es- 
sentially an infinity of false vacuum energy) and should 
have a horizon at some radius where 2m /r 1. Hence, 
it would seem that the critical behavior observed here is 
occurring within this horizon. 
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B. Addition of a regular scalar field 

An interesting result found by Choptuik in his study of 
collapse of a Yang-Mills field, is that the presence of even 
an infinitesimal amount of scalar field added to the Ein- 
stein Yang-Mills model drives the solution to the scalar 
field's critical solution at criticality Thus, the Yang- 
Mills field is unstable to the growth of the scalar field at 
criticality. 

To test for a similar effect here, I allow 'ip to be non- 
vanishing and look for the critical solution. Specifically, 
allowing the fields /(r, t) and '4){r, t) to be arbitrary Gaus- 
sian pulses, I can "mix" the amplitudes with a parameter 
a such that for cr = 0, the field if) vanishes and for cr = 1, 
the field / vanishes. 
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FIG. 4. Transition of critical solution from A = 0.46 to 
A = 3.44. A series of near-critical solutions are shown for 
a = 0.4 (solid), a = 0.5 (dotted), and a — 0.6 (dashed). 
The critical solution for small a is the new DSS, and as a 
increases, the critical solution becomes the original DSS. In 
the transition region, near-critical solutions appear to demon- 
strate structure from both critical solutions. 

For the EYM case, any initial data for which a ^ Q 
results in the scalar field critical solution being found . 
Here, no such sensitivity is observed. Instead, within 
some finite range of a the critical solution is observed to 
transition between the DSS of the scalar field (A = 3.44) 
and the new DSS presented here (A — 0.46). 

As cr is varied across the transition from the original 
DSS to the new one, the fields / and V' exchange dom- 
inance, both echoing at their respective value of A. In 
the transition region, the near-critical solution, as rep- 
resented by the field m'(r, i), displays this dual echoing 



and hence has a non-trivial structure. This transition is 
shown in Fig. ^. 

C. Critical Phenomena Inside a Monopole 

In contrast to the cases presented above in which 
/(oo,t) = 0, I consider initial data which have unit 
monopole charge as indicated by f{r,t) asymptoting to 
■q. Here, collapse within a monopole is modeled by first 
solving for the static monopole and then adding an in- 
going Gaussian perturbation to the profile /(r, t). This 
solution has the benefit of not being dynamic near the 
outer boundary. The Gaussian pulse is then parameter- 
ized by an amplitude p, and a search is conducted in p 
for a near critical solution. 

As stated, these evolutions should have the same criti- 
cal solution as those lacking the potential. This is verified 
numerically for rj € [0,0.2]. For rj large {rj > 0.20), no 
horizon-less, static solutions exist Q and hence the code 
is able to model these monopoles only within their hori- 
zons. A near-critical solution for rj = 0.15 is shown in 
Fig. I 
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FIG. 5. Near-critical solution obtained within a monopole 
for ri = 0.15. Notice that the field / asymptotes to ri for large 
r indicating the unit charge of the configuration. The discrete 
nature with A = 0.46 is observed to be the same as in the 
vanishing potential case. In comparison to the r] = case (see 
Fig. |l|), in the scaling regime (small r) the echoing behavior is 
evident, and in the large r regime the solutions are different. 

That static monopoles exist in this model naturally 
leads one to question whether any Type I critical phe- 
nomena exists. However, evolutions appear to indicate 
that the static monopoles with no horizon are not unsta- 
ble to collapse to a black hole. Hence, no such Type I 
critical phenomena has been observed. 
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IV. CONCLUSION 
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A new critical solution is found by assuming a hedge- 
hog ansatz for the triplet scalar field. The critical solu- 
tion exhibits discrete self-similarity with echoing period 
A = 0.46. The solution also exhibits power-law mass 
scaling with exponent 7 = 0.119. Results verify that the 
potential does not affect the critical behavior as expected. 

Furthermore, the threshold of black hole formation is 
examined within a global monopole. Though the po- 
tential is expected to play no role in Type II critical 
behavior, it does provide for the existence of static so- 
lutions, namely global monopoles. These static solutions 
then provide for the possibility that Type I behavior is 
present. The static solutions would have to be unstable 
to black hole formation for this Type I behavior to exist. 

However, evolutions do not show any instability of the 
monopole to collapse, and no Type I critical behavior is 
seen. A linear perturbation analysis should be sufficient 
to settle the question of their stability. If the static solu- 
tions are indeed stable, then no Type I behavior would 
be expected. 

Finally, as discussed in ||l^ symmetries of the ini- 
tial data can determine which critical solution is at- 
tracting. For the spherically symmetric complex scalar 
field, in addition to Choptuik's original critical solu- 
tion, another critical solution exhibiting a continuous 
self-similarity (CSS) exists An arbitrary family of 

initial data will find Choptuik's solution, not the CSS. 
However, certain types of initial data which maximize 
the [/(I) charge density are unable to shed their charge 
density and find instead the CSS. 

Here, instead of a complex scalar field (i.e. a doublet 
field) , the single scalar field is extended to a triplet field. 
Within the triplet scalar field model before making any 
assumptions on the triplet, at least two critical solutions 
exist, both DSS. The appropriate question is then will 
arbitrary families of initial data find Choptuik's original 
solution, the one presented here, or yet some other so- 
lution. To answer this with dynamical evolutions would 
require a general, three-dimensional collapse code. 

As a first step towards the goal of answering this ques- 
tion, Martin-Garcia and Gundlach ||l^ have shown that 
non-spherical perturbations to Choptuik's DSS decay 
and hence that at least some families of non-symmetric 
initial data will find Choptuik's DSS. However this work 
considered only a single scalar field, and does not address 
the case where a triplet is present. 

The question of the relative stability of the original 
DSS to this new DSS essentially reduces to determin- 
ing whether a family of generic, three-dimensional initial 
data for a triplet scalar field will shed its "in-phase" com- 
ponent to become the hedgehog DSS, or instead shed all 
its "out-of-phase" (i.e. non-zero monopole charge den- 
sity) component to become the scalar field DSS. 
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